# Metview Macro

t2m_forecast_grib = read("t2m_forecast.grib")
obs_bufr          = read("obs.bufr")


filter_obs_t2m = obsfilter(
    output : "geopoints",
    data   : obs_bufr
    )

fc_obs_diff = t2m_forecast_grib - filter_obs_t2m


shaded_coastlines = mcoast(
    map_coastline_land_shade        : "on",
    map_coastline_land_shade_colour : "cream",
    map_coastline_sea_shade         : "on",
    map_coastline_sea_shade_colour  : "HSL(193,0.78,0.94)"
    )

diff_symb_hot = msymb(
    legend                                 : "on",
    symbol_type                            : "marker",
    symbol_table_mode                      : "advanced",
    symbol_advanced_table_selection_type   : "list",
    symbol_advanced_table_level_list       : [0,1,5,10,15],
    symbol_advanced_table_max_level_colour : "red",
    symbol_advanced_table_min_level_colour : "white",
    symbol_advanced_table_colour_direction : "clockwise",
    symbol_advanced_table_marker_list      : 15,
    symbol_advanced_table_height_list      : [0.4,0.5,0.6],
    symbol_outline                         : "on",
    symbol_outline_colour                  : "charcoal"
    )


diff_symb_cold = msymb(
    legend                               : "on",
    symbol_type                          : "marker",
    symbol_table_mode                    : "advanced",
    symbol_advanced_table_selection_type : "list",
    symbol_advanced_table_level_list     : [-15,-10,-5,-1,0],
    symbol_advanced_table_colour_method  : "list",
    symbol_advanced_table_colour_list    : ["blue","sky","rgb(0.82,0.85,1)","white"],
    symbol_advanced_table_height_list    : [0.6,0.5,0.4],
    symbol_outline                       : "on",
    symbol_outline_colour                : "charcoal"
    )

plot(fc_obs_diff, diff_symb_hot, diff_symb_cold, shaded_coastlines)
